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Abstract 

At hadron colliders the dominant production mechanism of (be) mesons with large 
transverse momentum is due to parton fragmentation. We compute the rates and 
transverse momentum spectra for production of S-wave and P-wave (6c) mesons at the 
Tevatron via the direct fragmentation of the bottom antiquark as well as the Altarelli- 
Parisi induced gluon fragmentation. Since all the radially and orbitally excited (6c) 
mesons below the BD flavor threshold will cascade into the pseudoscalar ground state 
Be through electromagnetic and/or hadronic transitions, they all contribute to the 
inclusive production of Be. The contributions of the excited S-wave and P-wave states 
to the inclusive production of Be are 58 and 23%, respectively, and hence significant. 
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I. Introduction 


Bound states composed of two heavy quarks such as (cc), (bb) and (be) mesons, and (ccq), (bbq) 
and (beq) baryons (q denotes a light quark) are of interests because they can be produced 
with sizable rates at the current high energy hadron or e’'“e“ colliders. The Standard Model 
predictions of the production rates of these bound states can therefore be confronted with 
experimental data. 

The direct production of heavy mesons like heavy quarkonia and (fee) bound states can 
provide very interesting tests for perturbative QCD. The production of J/'0 and -0' [|^ at 
the Tevatron has already raised a lot of theoretical interests in explaining the excess of the 
experimental data above the lowest order perturbative QCD calculation |Q], especially at the 
large transverse momentum (px) region. The ideas of heavy quark and gluon fragmentation 
have been successfully applied to explain the experimental data of the prompt 
production from CDF within a factor of hve P JT(]|JTT|JT^ . Various attempts 
using the same fragmentation ideas have also been made to resolve the ip' surplus problem 
observed at CDF. 

Recently, the preliminary CDF results also showed that the production rates of the 
IS, 2S, and 3S T states are in excess of the leading order calculation. While the T(IS') 
and T(2S) results can be partly explained by including the fragmentation contribution, the 
T(3S') result showed an excess of about an order of magnitude over the QCD prediction 
even with the fragmentation contribution included |^. One subtlety is that the relevant pt 
for the fragmentation contribution to dominate should be larger in the bottomonium system 
than in the charmonium system, such that fragmentation is only valid for p^ > (1-2) nib in 
the former case. Besides, one has to worry about the very small pt region because, unlike the 
charmonium system, the experimental triggering conditions on the muon pair coming from 
the bottomonium leptonic decay allow experimentalists to measure the transverse momentum 
of the bottomonium all the way down to about 0-1 GeV [0. In order to fully understand 
the Pt spectrum of the bottomonium production, different production mechanisms have to 
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be brought into picture to explain the production rates in different pt regions. 

The {be) meson system, which is intermediate between the J/'ijj and T families, is also an 
interesting physical system to study. The mass spectrum of the (he) mesons can be predicted 
reliably from quarkonium potential models |T8|JT^ ,20[| without introducing any new param¬ 
eters and their decay constants can be computed using QCD spectrum sum rules [|I^pO| . 
We will adopt the Particle Data Group conventions, denoting the IS" pseudoscalar (^S'o) 
and vector meson (^S'l) (he) states by and S*, respectively. Higher radially and orbitally 
excited states are labeled by the standard spectroscopy notation: where the in¬ 

teger n is the principal quantum number, and L, S', and J = L + S are respectively the 
orbital angular momentum, total spin, and total angular momentum of the bound state. 
In the LS coupling scheme, for each principal quantum number n, the spin-singlet and the 
spin-triplet S-wave {L = 0) states are denoted by ^S'o and ^S*!, respectively. Those for the 
P-wave {L = 1) and D-wave {L = 2) states are denoted by ^Pi, ^Pj {J = 0,1,2), and ^D 2 , 
^Dj {J = 1,2,3), respectively. 

According to the results of the potential model calculation |]^ , the first two sets (n = 1 
and n = 2) of S-wave states, the first (n = 1) and probably the entire second set {n = 2) of P- 
wave states, and the first set (n = 1) of D-wave states lie below the BD flavor threshold. Since 
QCD interactions are diagonal in flavors, the annihilation channel of excited {be) mesons can 
only occur through the weak gauge boson (W) exchange and is therefore suppressed relative 
to the electromagnetic and hadronic transitions to other lower lying states. The excited 
states below the BD threshold will cascade down into the ground state Be via emission of 
photons and/or pions, while the other states above the BD threshold will decay rapidly into 
a pair of B and D mesons. Inclusive production of the Be meson therefore includes the 
production of the n = 1 and n = 2 S-wave and P-wave states, and the n = 1 D-wave states. 

The production of the S-wave and ^Si) states were first computed exactly to leading 
order in Ref. at the machine, in particular at the Z resonance. Later, it was 

realized |^,§] that the complicated formulas in these complete calculations can be simplified 
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by a factorization approach. The dominant contribution in the leading order calculation can 
be factorized into a short distance piece, which describes the partonic process of the decay 
of Z into a high energy bb pair, and a fragmentation piece describing how the b antiquark 
splits into the two S-wave states. The corresponding fragmentation functions and 

which are independent of the short-distance piece, were shown to be calculable 
by perturbative QCD at the heavy quark mass scale [Q. Recently, the production of the 
S-wave states has also been computed at hadronic colliders like the Tevatron and the Large 
Hadron Collider (LHC) both by a complete 0(ag) calculation [pSjpH 


and by using 


the simpler fragmentation approach 


, P^] . We note that unlike the J/ip production in 
hadronic collisions in which the major contributions come from (1) g —>■ Xcj fragmentation 
followed by the decay XcJ J/'ip + x [Pi TO P ^ P^, and (2) gluon fragmentation into a 


color-octet (cc) state which subsequently evolves nonperturbatively into J/V' [pRtpR 


the fragmentation diagrams for producing (be) mesons form a subset of the whole set of 
0{ag) diagrams, which are the leading-order diagisons for producing (fee) mesons in hadronic 
collisions. It is therefore not clear that if the production of (fee) mesons in hadronic collisions 
is dominated by parton fragmentation. However, detailed calculations by Chang et al. [P3| , 
Slabospitsky [^, and Kolodziej et al [^] indicated that the fragmentation approach is valid 
for the S-wave production at the large transverse momentum region. We will discuss more 
about this later in the closing section. 

In this paper, we study the hadronic production of (fee) mesons. We compute the pro¬ 
duction rates of the S-wave and P-wave (fee) mesons at the Tevatron using the fragmentation 
approach. Intuitively, the dominant production mechanism of the (fee) mesons at the large 
transverse momentum region must arise from the direct fragmentation of the heavy fe anti¬ 
quark. The relevant question is whether experiments can probe the transverse momentum 
region where fragmentation dominates. Unfortunately, to answer this question it also re¬ 
quires a complete 0{a'^) calculation for the production of the P-wave states, which is not 
available at the moment. Here we assume that the fragmentation approach also works for 
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production of the P-wave {be) mesons in the transverse momentum range that is being 
probed experimentally at the Tevatron. In this work, we do not include the contributions 
from the D-wave states because the corresponding fragmentation functions are very small 
E3. Although the production of the D-wave states is of great interest by themselves, they 


only contribute about 2 % |^| to the inclusive production of Be- In the near future, like 
other heavy quarkonia, the production of {be) mesons at the Tevatron may therefore provide 
another interesting test for perturbative QCD. Although we will only show our results for 
the positively charged states {be), all the results presented in this work also hold for the 
negatively charged states {be). 

The organization of this paper is as follows. In the next section, we discuss in detail 
the general procedures to calculate the production cross sections using the fragmentation 
approach. In Sec. Ill we present the transverse momentum spectra and the integrated 
cross sections for the production of S-wave and P-wave {be) states, as well as the inclusive 
production rate for the Be meson. Discussions and conclusions are made in Sec. IV. For 
completeness we also collect all the S-wave and P-wave fragmentation functions at the heavy- 
quark mass scale in the Appendix. 

Before leaving this preamble, we note that a preliminary result from CDF had provided 
a hint for the Be existence by looking at the production rate of J/ip + tt in the mass bin of 
6.1-6.4 GeV M- 


II. Inclusive production cross sections in the fragmen¬ 
tation approach 

Theoretical calculations of production cross sections in high energy hadronic collisions are 
based on the idea of factorization. Factorization divides an inclusive or exclusive hadronic 
production process into short-distance pieces and long-distance pieces. The short-distance 
pieces are perturbatively calculable to any desired accuracy in QCD, while the long-distance 
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pieces are in general not calculable within perturbation theory but can be parameterized 
as phenomenological functions, which can be determined by experiments. The factorization 
used here for the production of (be) mesons divides the process into the production of a 
high energy parton (a b antiquark or a gluon) and the fragmentation of this parton into 
various (be) states. The novel feature in our approach, which is due to a recent theoretical 
development |^,^], is that the relevant fragmentation functions at the heavy quark mass scale 
can be calculated in perturbative QCD to any desired accuracy. This is easily understood 
from the fact that the fragmentation of a 6 antiquark into a (be) meson involves the creation 
of a cc pair out of the vacuum. The natural scale for this particular hadronization is of order 
of the mass of the quark pair being created. In our case, this scale is of order rric, which 
is considerably larger than Aqcd- One can therefore calculate reliably the fragmentation 
function as an series expansion in the strong coupling constant using perturbative QCD. 

The production of high energy partons also involves the factorization into the parton 
distribution functions inside the hadrons and the parton-parton hard scattering. Let H 
denotes any (be) meson states. The differential cross section da/dpT versus the transverse 
momentum pt of H is given by 


da 

dpT 


(pP^H(pt)X) = ^ / (ixida;2(i^/i/p(a;i,/i)/j/p(a;2,/i) 




dd 

dpT 


(ij b(pT/z)X, p) 


dd 


xDi^h(Z:P) + ^ 9(Pt/z)X,p) Dg^H(z,p) 


( 1 ) 


The physical interpretation is as follows: a heavy b antiquark or a gluon is produced in a 
hard process with a transverse momentum Pt/z and then it fragments into H carrying a 
longitudinal momentum fraction z. We assume that H is moving in the same direction as 
the fragmenting parton. In the above equation, /j/p(p)(a;,/x)’s are the parton distribution 
functions, dd’s represent the subprocess cross sections, and Di^H(z, pYs are the parton 
fragmentation functions at the scale p. For production of the b antiquark, we include the 
subprocesses gg —>• bb, gb —>• gb, and qq —>• bb; while for the production of the gluon g, we 
include the subprocesses gg —>• gg, qq —> gg, and gq(q) —>■ gq(q)- In Eq. (0), the factorization 
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scale /i occurs in the parton distribution functions, the subprocess cross sections, and the 
fragmentation functions. In general, we can choose three different scales for these three 
entities. For simplicity and ease of estimating the uncertainties due to changes in scale, we 
choose a common scale /i for all three of them. 

The physical production rates should be independent of choices of the scale /i, because fi 
is just an artihcial entity introduced to factorize the whole process into different parts in the 
renormalization procedure. However, this independence of scale can only be achieved if both 
the production of the high energy partons and the fragmentation functions are calculated to 
all orders in as- So far, only the next-to-leading order d’s and the leading order fragmentation 
functions are available, so the production cross sections do depend on the choice of /x to a 
certain degree. We will estimate the dependence on /x by varying the scale /i = (0.5 — 
where fiR is our primary choice of scale 

Hr = (parton) + mf . (2) 

This choice of scale, which is of order p'p(parton), avoids the large logarithms in the short- 
distance part (j’s. However, we have to sum up the logarithms of order HR/^b in the frag¬ 
mentation functions. But this can be implemented by evolving the Altarelli-Parisi equations 
for the fragmentation functions. 

11, a Evolution of Fragmentation Functions 

The Altarelli-Parisi evolution equations for the fragmentation functions are 

^ rl rl 

—Pb^b{z/y,H)Di^H{y,H)+ —Pb^g{z/y,H) Dg^H{y,H), (3) 

ofi Jz y Jzy 

d dy dy 

H^Dg^H{z,H)= —Pg^i{z/y,H)Di^H{y,fx)+ —Pg^g{z/y,H)Dg^H{y,H), (4) 

oy, Jzy Jzy 

where H denotes any {he) states, and Pi^j are the usual Altarelli-Parisi splitting functions. 
The leading order expressions for Pi^j can be found in the Appendix. 
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The boundary conditions for solving the above Altarelli-Parisi equations are the fragmen¬ 
tation functions ^q) and Dg^^iz, ^q) that we can calculate by perturbative QCD 

at the initial scale fiQ, which is of the order of the &-quark mass. At present, all the S-wave 
0 and P-wave [0 fragmentation functions for b (be) have been calculated to leading order 
in as- They are all collected in the Appendix. The initial scale po for Di^jj{z, is chosen 
to be po = '^b + 2mc, which is the minimum virtuality of the fragmenting b antiquark [Q. 
On the other hand, since the initial gluon fragmentation function Dg_,H{z, po) is suppressed 
by one extra power of as relative to fio), we simply choose the initial gluon frag¬ 

mentation function to be Dg^niz, fJ-o) = 0 for po < 2(mb -|- rric) — the minimum virtuality 
of the fragmenting gluon [^. We can also examine the relative importance of these frag¬ 
mentation functions. The initial fio) is of order whereas the initial Dg^niz, fio) 

is of order af and has been set to be zero for po < 2(mb -|- rric) as discussed above. But 
when the scale is evolved up to a higher scale /i, Di^jj{z, fi) is still of order while the 
induced Dg^H(^z, fi) is of order a^log(p/po)- At a sufficiently large scale fi the logarithmic 
enhancement can offset the extra suppression factor of as- Thus the Altarelli-Parisi induced 
gluon fragmentation functions can be as important as the b antiquark fragmentation, even 
though the initial gluon fragmentation functions are suppressed. While these Altarelli-Parisi 
induced gluon fragmentation functions play only a moderate role at the Tevatron, they will 
play a more significant role at the LHC PU| . 

To obtain the fragmentation functions at an arbitrary scale greater than fiQ, we numeri¬ 
cally integrate the Altarelli-Parisi evolution equations (|^)-(^ with the boundary conditions 
described in the above paragraph. Since the initial light-quark fragmentation functions 
Dq^ni^z, fio) are of order one can set them to be zero as well for po < 2(mb -|- rric)- 
One may ask if the light-quark fragmentation functions can be induced in the same manner 
as in the gluon case by Altarelli-Parisi evolution. Since both the Altarelli-Parisi splitting 
functions Pq^q and Pq^g are total plus-functions, the induced light-quark fragmentation 
functions Dq^niz) can only be total plus-functions or vanishing identically. We do not an- 


ticipate that these induced light-quark fragmentation functions, if nontrivial, will play any 
signihcant role in our analysis. 


III. Numerical results 


Leading order QCD formulas are employed for the parton-level scattering cross sections 
and CTEQ(2M) is used for the initial parton distributions. The inputs to the initial 
fragmentation functions are the heavy quark masses rrib and rUc, and the nonperturbative 
parameters associated with the wavefunctions of the bound states. For the S-wave states 
there is only one nonperturbative parameter, which is the radial wavefunction i?n 5 ( 0 ) at 
the origin. However for the P-wave fragmentation functions we have two nonperturbative 
parameters Hi and H'^ associated with the color-singlet and the color-octet mechanisms, 
respectively |^. Two of the P-wave states (^Pi and ^Pi) are mixed to form the two physical 
states, denoted by and ll’*"'). Further details of the mixings can be found in Refs. [[7tp!8|. 
The above input parameters for the fragmentation functions are summarized in Table 

As mentioned in the previous section, we have set the scale p in the parton distribution 
functions, subprocess cross sections, and the fragmentation functions to be the same. We will 
later vary /i between 0.5/1^ and 2/iij, where hr is given in Fq. (^), to study the dependence 
on the choice of scale. For the strong coupling constant as(p) entered in the subprocess cross 
sections, we employ the following simple expression 




as{Mz 


(5) 


1+ H^«.(M.)ln(^)’ 
where u/ is the number of active flavors at the scale /i and as{Mz) = 0.118. In order to 
simulate the detector coverage at the Tevatron, we impose the following acceptance cuts on 
the transverse momentum and rapidity of the {he) state H : 


PT{H)>6GeV and \y{H)\ < 1 . 


( 6 ) 


The numerical results of the pr spectra for the {be) state H with various spin-orbital 
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quantum numbers are shown in Fig. and Fig. ^for the cases of principal quantum number 
n = 1 and n = 2, respectively. The integrated cross sections versus are also shown 

in Figs. ^ and ^ for the cases of n = 1 and n = 2, respectively. 

Now we can predict the inclusive production rate of the meson. As the annihilation 
channel is suppressed relative to the electromagnetic/hadronic transitions, all the {be) excited 
states below the BD threshold will decay eventually into the ground state B^ via emission 
of photons or pions. Since the energies of these emitted photons and/or pions are limited by 
the small mass differences between the initial and hnal (he) states, the transverse momenta 
of the (he) mesons will not be altered appreciably during the cascades. Therefore, we can 
simply add up the pt spectra (Fig. ^ and Fig. |^) of all the states to represent the pt spectrum 
of the inclusive Be production. Similarly, we can add up the integrated cross sections (Fig. ^ 
and Fig. ^ to represent the integrated cross section for the inclusive Be production. Thus, 
we can obtain the inclusive production rate of Be as a function of p^™{Be). Table gives the 
inclusive cross sections for the Be meson at the Tevatron as a function of p™™(i?c), including 
all the contributions from the n = 1 and n = 2 S-wave and P-wave states. These cross 
sections should almost represent the inclusive production of Be by fragmentation, because 
the contributions from the D-wave states are expected to be minuscule. The results for 
h = hr?/2, Pr, and 2pii are also shown. The variation of the integrated cross sections with 
the scale p is always within a factor of two, and only about 20% for p-r > 10 GeV. We will 
discuss more about the dependence of scale in the next section. 

At the end of Run Ib at the Tevatron, the total accumulated luminosities can be up 
to 100-150 pb“^ or more. With px > Q GeV, there are about 5 x 10® B(' mesons. The 
lifetime of the Be meson has been estimated to be of order 1-2 picosecond |^, which is 
long enough to leave a displaced vertex in a silicon vertex detector. Besides, Be decays into 
J/ih + X very often, where X can be a tt’*', p"*", or and J/(j can be detected easily 
through its leptonic decay modes. The inclusive branching ratio oi Be ^ J/ip + X is about 
10% IQ. When X is e+z/g or we will obtain the striking signature of three-charged 
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leptons coming off from a common secondary vertex. The combined branching ratio of 
Be —>■ {i,i' = e,/i) is about 0.2%. This implies that there will be of 

order 10^ such distinct events for 100 pb“^ luminosity at the Tevatron. Even taking into 
account the imperfect detection efficiencies, there should be enough events for conhrmation. 
However, this mode does not afford the full reconstruction of the B^. If X is some hadronic 
states, e.g., pions, the events can be fully reconstructed and the B^ meson mass can be 
measured. The process > J/'ijj + 7 t^ —> is likely to be the discovery mode for B^. 

Its combined branching ratio is about 0.03%, which implies about 300 such distinct events 
at the Tevatron with a luminosity of 100 pb“^. 

After the next hxed target runs at the Tevatron, the Main Injector will be installed in 
1996-1997 according to the present plan [^]. The Main Injector will give a signiheant boost 
in the luminosity while the center-of-mass energy stays the same. The upgraded luminosity 
is estimated to be about ten times larger than its present value. This enables Run II to 
accumulate a total luminosity of 1-2 fb“^, which implies that about 10^ —10® Be mesons will 
be produced. With the Main Injector installed, it might be possible to produce the D-wave 
(be) states with sizable rates. 

Tevatron will continue running until the next generation of hadronic colliders, e.g., the 
LHC. The present design of the LHC is at a center-of-mass energy of 14 TeV and the yearly 
luminosity is of order 100 fb“^. In Table 13 we show the inclusive cross sections for the Be 
meson at the LHC as a function of including the contributions from n = 1 and n = 2 

S-wave and P-wave states. With the assumed 100 fb”^ luminosity there are about 3 x 10® 
Be mesons with > 10 GeV and \y(Be)\ < 2.5. With such a high luminosity at LHC, one 
expects sizable numbers of the various D-wave (be) states to be produced as well. The LHC 
will then be a copious source of (be) mesons such that their properties e.g., spectroscopy 
and decays, can be thoroughly studied. In addition, the mixing and CP violation studies 
are possible. For example, one can use Be to tag the flavor of Bg in the decay B^ —> B^i~^i/i 
for the studies of R® — R® mixing. Also, CP violations in the Be system can be studied by 
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looking at the difference in the partial decay widths of B'^ X and X. 

IV. Discussions and conclusions 

First, we briefly discuss the various sources of uncertainties in our calculation. One uncer¬ 
tainty conies from the use of the naive Altarelli-Parisi evolution equations. As pointed out 
, the naive Altarelli-Parisi equations in Eqs. (^-(||) do not respect the phase 
space constraints. Inhomogeneous evolution equations were then advocated to remedy for 
these problems. The major effect is to correct the unphysical blow-up of the evolved gluon 
fragmentation functions when z gets too close to 0. The corrected gluon fragmentation func¬ 
tions, instead of having the unphysical blow-up at z = 0, turn to zero smoothly below a 
certain threshold value of z. Despite the dramatic changes of the evolved gluon fragmenta¬ 
tion functions for small z values, this effect will not show up easily in our calculation because 
the small z region is very likely excluded by the transverse momentum cut imposed on the 
(be) mesons. Therefore, in this paper we keep on using the homogeneous Altarelli-Parisi 
equations in Eqs. (3)-® to evolve our fragmentation functions. 

A second source of uncertainty is due to the choice of the factorization scale p. We show in 
Fig. ^the dependence of the differential cross sections on the factorization scale /i by plotting 
the results for the various choices of /i = and 2/rjj. For clarity we only show the 

curves for the 1 ^S'o and 1 states in Fig. |^. The behaviors for other states are similar. Note 
that in the /i = Hr/ 2, fiR, and 2^r curves, the running scales used in the strong coupling 
constant which entered in dd’s and in the parton distribution functions fi{x)^s are equal 
to /i, while the running scale used in the fragmentation functions is set to be max (/i, (Iq), 
where po is the prescribed initial scale for the fragmentation functions. Figure ^ shows 
that the change in the factorization scale /i gives different results for the differential cross 
sections. The fi = 2 ^r curves show that the differential cross sections increase (decrease) 
only slightly at the low (high) pr region. Although the /i = /iij/2 curves demonstrate 


in Refs. 
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larger changes in the differential cross sections, the variations are always within a factor of 
two. The integrated cross sections at varions as already shown in Table y, also 

indicate lesser sensitivity in the scale as one increases The variations of differential 

cross sections with the scale /i demonstrate the effects of the next-to-leading order (NLO) 
corrections. Only if all the NLO corrections are calcnlated for each piece contained in Eq. 
namely the parton distribntion fnnctions, the parton-parton hard scattering cross sections, 
and the fragmentation fnnctions, can these variations be rednced snbstantially. Since the 
pertnrbative QCD fragmentation fnnctions are only calculated to leading order, the NLO 
calculations of the perturbative QCD fragmentation functions can provide an improvement to 
our calculation. However, due to the rather weak dependence on the scale in our leading order 
calculation, our results should be rather stable under higher order perturbative corrections. 

Other uncertainties come from the input parameters to the boundary conditions of the 
fragmentation functions. These are the heavy quark masses mb and me, and the nonpertur- 
bative parameters describing the bound states. Slight changes in me and mt could possibly 
lead to appreciable changes in the fragmentation functions, as indicated by the m^ and 

dependence, respectively in the denominators of the S-wave and P-wave fragmentation 
functions. (Note that Hi/m ~ 9|i?'(0)p/(327rm^).) However, in the numerators the wave- 
functions at the origin |i?(0)p and |i?'(0)p also scale like m^ and m®, respectively. Therefore 
the dependence on the heavy quark masses should be mild in the fragmentation functions. 
The color-octet parameters HgS are associated with the additional color-octet contributions 
and they are not well determined. However, we do not expect that they will play any signif¬ 
icant role in the present context and only refer our readers to Ref. [0] where some discussions 
of these parameters can be found. 

There was a controversy in the recent literature concerning about the 

importance of parton fragmentation in the production of the Be {B*) meson at hadronic 
colliders. The controversy arises because the Feynman diagrams responsible for the b anti¬ 
quark fragmentation form a subset of the whole set of Feynman diagrams, which contribute 
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at the order of Thus, there is a competition between the fragmentation contribution and 
the non-fragmentation contribution. Some authors had referred to the latter contribution 

have presented such a 


as recombination. So far, hve independent groups 
complete 0{aj) calculation. Chang et al. and Slabospitsky ||^ agreed that the frag¬ 
mentation contribution dominates at the large transverse momentum region. However, we 
could not hnd in their work the precise value of pt at which the fragmentation contribution 


begins to dominate. On the other hand, Berezhnoy et al. [^| claimed that fragmentation 
never dominates for all kinematical region of pt and the recombination diagrams can never 


be ignored. Independently, Masetti and Sartogo |26] have recently obtained results claimed 


to be consistent with Berezhnoy et a/ [^|. There are certainly discrepancies among these cal¬ 
culations. Most recently, Kolodziej et al performed yet another independent calculations 
using two different methods - the usual trace and the modern helicity amplitude techniques. 
Both methods agreed with each other numerically to very high accuracy, according to the 
authors. Kolodziej et al. also did a very thorough comparison with previous exact cal¬ 
culations by comparing numerically their matrix elements with other groups P3| , |2^P3[] , but 
could not hnd agreement using identical input parameters. Furthermore, Kolodziej et al. 
P7| did a comparison of the production of Be mesons between the exact calculation and the 
fragmentation approximation, and found that the fragmentation approach is indeed valid for 
Pt > 10 GeV! Because of the conclusions of Ref. [^, we believe that the fragmentation con¬ 
tribution should begin to dominate the production of (he) mesons at px ^ 10 GeV. Also, we 
would like to emphasize that the detector coverage and performance do require a minimum 
transverse momentum cut and a rapidity cut on the Be meson. Thus, the non-fragmentation 
(recombination) contribution, which is substantial at low px, becomes less important as the 
lowest px range is being excluded. 

Tn their most recently revised versions, both groups of Chang et al. (the second entry of Ref. Q) and 
Berezhnoy et al. confirmed the results of Kolodziej et al. . Thanks to an anonymous referee for 
pointing out the situation. 
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The main result of this paper is the inclusive cross section of the ground state B^. at the 
Tevatron. We have included in this work all the contributions from the direct b fragmentation 
and induced gluon fragmentation, as well as the cascade contributions from all the excited 
S-wave and P-wave states below the BD threshold. Previous result presented by one of 
us in Ref. 1 28[| only included the direct b fragmentation into S-wave states, since the P- 
wave fragmentation functions were not available at that time. Later it was pointed out in 
Ref.[p9| that the induced gluon fragmentation contribution is nonnegligible. The induced 
gluon fragmentation contribution is about 20 - 30% (30 - 40%) of the direct b fragmentation 
contribution at the Tevatron (LHC) energies. The two sets of P-wave states included in this 
work contribute at the level of 30%. Furthermore, we have employed here a more updated 
set of parton distribution functions than those used in Refs.. As a result, the cross 
sections increase roughly by 25%, mainly due to the increase in the gluon luminosities. When 
all these aforementioned differences are taken into account, the net result is about 70 - 80% 
(a factor of two) larger than the previous result presented in Ref. PH| at the Tevatron (LHC). 

To recap, in this paper we have performed calculations of the transverse momentum 
spectra and integrated cross sections for the S- and P-wave (be) mesons. We have also 
predicted the inclusive production rate for the Be meson, and found that there should be 
enough signature events to conhrm the existence of Be at the Tevatron, whereas the LHC 
will be a copious source of Be- Nevertheless, one should keep in mind that what we have 
calculated in this paper represents only the contribution from fragmentation, while there 
should also be other non-fragmentation contributions, especially at the low pt region (pt ^ 
10 GeV). At least, our results represent a lower bound of the production rate for Be. In the 
near future, the production of Be at the Tevatron would provide another interesting test for 
perturbative QCD, and we expect a very exciting experimental program of Be at the LHC. 
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Appendix 


The Altarelli-Parisi evolution equations were given in Eqs. 
list all the Altarelli-Parisi kernels in the following, 

4as(/i) f 1 + x'^\ 


Pl^g{x,fi) = 
Pg^i{X,fi) = 


Svr \1 — X J ^ ' 
4Q;^(/i) / 1 + (1 - xf ' 

Svr y X 
as{^) 

271 

6as(/i) 


+ (1 - x)2) , 

1 — X 


X 


+ 


TT 


X 


). For convenience, we also 


(A.l) 

(A.2) 

(A.3) 

(A.4) 


+ x{l - x) + [^ - ^] 6{1 - x) 


(1 -x)_ 

where uj is the number of active flavors below the scale /r. 

The fragmentation functions fio) for the S-wave and P-wave (be) states have 

been obtained to leading order in strong coupling constant a* at the heavy quark mass scale 
Po = + 2mc- We introduce the following notations: r = mc/{mb + me), f = 1 — r, and 

m = mbmc/(mb + me). The nonperturbative parameters are i?„s(0) for the S-wave states, 
and Hi(n) and H'^(n) for the P-wave states. The mixing angle between the ^Pi and ^Pi states 
is denoted by cos 9 rr,w . The numerical values of these input parameters can be inferred from 
potential model calculations (see for example |T^) and are listed in Table | for convenience. 
The scale of the strong coupling constant entered in Di^u(z, ^q) is set to be 2mc- 
The S-wave fragmentation functions are ||^ 

2a^(2mc)^|i?ns(0)P rz(l-z)‘^ 


>6c(nlSo) 


X 


SlTTm^ (1 — rz)^ 

6 - 18(1 - 2r)z + (21 - 74r + mr‘^)z‘^ 


2r(6 — 19r -|- 18r^)z^ -f 3r^(l — 2r -|- 2r^)z 


4 


(A.5) 
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^b^bc{n^S\) /^o) 


2as{2mcY\Rns{^)\^ rz{l - zf 


X 


277rml (1 — rz)® 

2{3 — 2r)z + 3(3 — 2r + 4r^)2:^ — 2f(4 — r + 2r^)z^ 


+ r^{3 — 2r + 2r'^)z^ 


(A.6) 


Polarized S-wave fragmentation functions can also be found in Ref. but they are not 
used in this work. 

The P-wave fragmentation functions consist of two terms, the color-singlet piece {Hi 
term) and the color octet piece (ifg term) [^: 


^b—>bc{n^Pl) h'o) 
^b—^bc{n^Pj) {^1 h'o) 


Hi(bc)in) 


m 

Hi(bc)M 


fe^6c(iPi)(^’ ^) + 3 


^8(5c)(^)^(8) 


m 


fe^fec(iSo) 


» , 


(A.7) 


-D 


b^bc{3Pj)i^^ + (2J + 1) 


m 


-^8(5c)(^) ^(8) 


m 


fe^5c(3Si) 


», 


(A.8) 


where A is a factorization scale for the P-wave fragmentation functions. It separates the 
reduced heavy quark mass scale m and the scale mv which characterizes the bound state 
structures with v being the typical relative velocity of the heavy quark and antiquark inside 
the bound state. In our numerical work, we will choose A ~ m. 

At leading order the color-singlet pieces do not depend on A and they are given by 


-Dfe2.6c(iPi)(^) 


16Q!^(2mc) rr^z{l — zY 


6 — 6(4r^ — 8r -|- 5)z 


243 (1 - fz) 

+ (32r^ - 96r^ + 250r2 - 210r + 69)^^ 

-1- 8f(4r^ -I- 12r^ — 48r^ -|- 37r — 12)^^ 

+ 2f2(16r^ 4- leir^ - 114r + 42)/ - 6/(4/ + 28/ - 15r + 7)/ 
-1- f‘^{4Qr‘^ — l4r + 9)/ 


(A.9) 
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D 


W_ 

b^bc{^ Pq) 


(^) 


16a‘^{2mc) rr^z{l — z) 


6(4r-l)2 + 6(4r-l)(20r2 


729 (1 - fz)^ 

(832p - 1456r^ + lOSSr^ - 362r + 63);^^ 


16r + 5);^ 


- 8f(100r^ - 184r^ + 118r2 - 22r + 9);^^ 

+ 2f2(416r^ - 776r^ + 369r2 + 42r + 24)/ 

- 2/(240/ - SlOr^ + 232r2 + 59r + 9)/ 
+ /(96/- 240/+ 134/+ 34r + 3)/ , 


(A.IO) 


L 


= 


32Q;/2mc) rr^z(l — z) 


and 


6 + 6(4r — 5)z 


243 (l-f/8 

+ (lOr"^ + 64r^ — 98r + 63)/ + 8f{2r* + 2r^ — 13r^ + Hr — 9)/ 

+ 2/(8r"^ — 16r^ + 47r^ — 18r + 24)/ + 2/(8r^ — 24r^ + r — 9)^:^ 

+ /(12r2 + 2r + 3)/ , (A.ll) 


I. 


b^bc(^P2) 


(/ = 


64Q;/2mc) rr^z{l — z')^" 


12 + 12(2r - 5);^ 


729 (1 - fzf 

+ (92r2 - 76r + 135)/ + 4(10r3 - 54r2 + 31r - 45)/ 

+ 2(46r^ - 16r=^ + 123r2 - 78r + 75)/ 

— 4f(6r"^ + 9r^ + 40r^ — 13r + 18)/ 

+ /(12r^ - 12r=^ + 55r2 - lOr + 15)/ ' 


(A. 12 ) 


The ^Pi and ^Pi states are mixed in general, giving rise to the following physical mass 
eigenstates |/') and ll’*'), 


ll"') 

= cosOmipPi) - sin 6 'mix|^Pi) , 

(A. 13 ) 

|V) 

= sin^mix^A) + cos 6 'j„ix|^Pi) . 

(A. 14 ) 
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Thus, in general, we have 






cos' (^) + sin' 

sin 0jnix COS 0mix-Djnix(^) ) 
sin^ 0mixT>^2bc(lPi)(^) + COS^ ^mix^^2fe,(3p,)(^) 
T sin cos , 


i(i) 


b^bc(l+ ) 


o™Sc,i*)(‘) 


(A.15) 


(A.16) 


with 


(1) 323/2a^(2mc) rr^z(l - 2:)" 

-^mixl^J ~ 030 /I 


6 — 6(2r^ — 4r + 3);^ 


243 (1 - fz)^ 

(24r^ + 52r^ — 52r + 21)2:^ + 2f(14r^ — 6r^ + 15r — 6)2:^ 

f^(2r^ + 8r — 3)2:"^"^ 


(A.17) 


Finally, the color-octet pieces are given by 




a^(2mc) rr^z(l — z)^ 


6 - 18(1 - 2r)z + (21 - 74r + 68r^)z^ 


162 (1 - fz)^ 

— 2f(6 — 19r -|- 18r^)2:^ -f 3f^(l — 2r -|- 2r^)z^ 


(A.18) 


and 


b£so,.s)(") 


a^(2mc) rr^z(l — z) 


2 - 2(3 - 2r)z + 3(3 - 2r + 4r^)z^ 


162 (1 - fz)^ 

— 2f(4 — r -T 2r^)z^ + f^(3 — 2r -|- 2r^)z^ 


(A.19) 
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Table I: Input parameters to the perturbative QCD fragmentation functions for n = 1 and n = 2. 



n = 1 

n = 2 

rUb 

4.9 GeV 

4.9 GeV 

rric 

1.5 GeV 

1.5 GeV 

-Rns(O) 

1.28 GeV3/2 

0.99 GeV3/2 

Hi 

10 MeV 

14 MeV 

H'M) 

1.3 MeV 0 

1.8 MeV 0 

COS ^mix 

0.999 

0.957 


Table II; The inclusive production cross sections for the meson at the Tevatron including the 
contributions from all the S-wave and P-wave states below the BD threshold as a function of 
The acceptance cuts are pt{Bc) > 6 GeV and \y{Bc)\ < 1. 


(GeV) 


a (nb) 




p = Pr 

fJ' — 2/ii? 

6 

2.81 

5.43 

6.93 

10 

0.87 

1.16 

1.22 

15 

0.26 

0.29 

0.26 

20 

0.098 

0.097 

0.083 

30 

0.021 

0.018 

0.014 
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Table III; The inclusive production cross sections for the meson at the LHC including the 
contributions from all the S-wave and P-wave states below the BD threshold as a function 
of The acceptance cuts are pt{Bc) > 10 GeV and \y{Bc)\ < 2.5. 


p?l“ (GeV) 

a (nb), p = Pr 

10 

33.8 

20 

4.4 

30 

1.1 

40 

0.41 


Figure Captions 

1. The differential cross section da/pT versus pr of the {be) meson {H) in various spin- 
orbital states with n = 1 at the Tevatron. The acceptance cuts are Pt{H) > 6 GeV 
and \y{H)\ < 1. 

2. Same as Fig.l for n = 2. 

3. The integrated cross section a{pT > p™™) versus the minimum p™™ cut on the {be) 
meson {H) in various spin-orbital states with n = 1 at the Tevatron. The acceptance 
cuts are the same as Fig.l. 

4. Same as Fig.3 for n = 2. 

5. The comparison of the differential cross sections da/dpr at different factorization scales 
for the 1 ^5*0 and 1 states at the Tevatron. The acceptance cuts are the same as 

Fig.l. 
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